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ABSTRACT 


Different  forms  of  correlated  atomic  wave  func¬ 
tions  and  their  implications  on  the  integrals  occurring 
in  calculations  of  variational  solution  to  Schr&dinger' s 
equation  are  briefly  discussed*  A  scheme  for  computing 
the  integral 

is  described  and  applied  to  the  case 


I.  INTRODUCTION 


In  ntoale  and  anlneular  thaory  aavaral  variational  tolu> 
tlona  to  tha  Sehroadlngar  aquation  hava  baan  abtalnad  for  a  varl- 
aty  of  nMll  ayataaa  ualng  trial  wava  funetlona  of  a  varying  degree 
of  coaplaxlty.  Tha  dlffarant  Mthoda  of  calculating  variational 
■olutlons  ara  coaawnly  claaatflad  according  to  tha  analytical  form 
of  trial  wava  function  uaad.  Tha  mathod  of  so>callad  correlated 
wava  functiona.  i.a.  funetlona  that  dayand  aspliclt|y  on  tha  Inter- 
alactronlc  dlataneaa  Cl  * |  Ij  •  ^  I  ***  first  Introducad  by 
Hyllaraaa.  Tha  method  has  latar  bean  usad  in  aumsrous  applications 
on  small  systams,  pra'doailnantly  two-alactron  systems.  Ona  of  the 


1)  P.O.  L<hfdln»  Adv.  Cham  Phys.  2  ,  207  (Intarsclanca  Publ.  L 

Prigogina.Ed.  New  York,  1959) 

2)  B.A.  Hylleraas,  Z.  Phys.  ^  ,  347  (1929) 


amrits  of  the  amthod  Is  that  it  gives  good  results  and  yat  can 

give  wava  functions  in  an  attractively  condensed  fora.  It  seems 

to  be  fairly  difficult,  however,  to  ganaraliae  the  method  to 

systems  with  sure  than  two  electrons.  Ona  of  the  main  obstacles  is 

that  considerabla  difficulties  are  mat  in  calculating  tha  integrals 

with  several  interalectronic  distances  in  tha  integrand.  Recently 

soma  progress  has  h*an  made  in  generalising  tha  Mthod  of  correlated 

3) 

wava  functions  to  many  electron  atoois,  but  still  an  accurate  calcu> 
lation  of  the  integrals  containing  functiona  of  seams  to  he 
an  unsolved  problem  in  its  full  generality. 


Savaral  diffarant  Mthoda  hava  baan  auggaatad  how  to 


3)  L.  Siaaa.  Fhya.  Rav.  m  ,  169  (1962) 

4)  H.M.  JaMa  and  A.S.  Coolidga,  Phya.  Rav.  ^9  688  (1936) 

5)  V.  Pock,  H.  Vaaaalov  and  M.  Patraahan,  J.  Bxplo.  Theorat.  Phya. 

(USSR),  ^  ,  723  (1940) 

6^tP«  Walah  and  8.  larawita,  Phya.  Rav.  m  ,  1206  (1959) 

7)  L.  Saaaa,  J.C.P.  »  ,  1072,  (1961) 

8)  J.L.  Calala  and  P.O.  Ltfwdln  ,  J.  Mol.  Spac.  8,  203  (1962) 


calculata  ouch  intagrala  and  tha  purpoaa  of  thia  note  la  to 
daacriba  a  alapla  achaaw  for  coag>utlng  cartain  atoadc  Intagrala 
that  hava  an  axpliciC  dependence  of  lntar>alaetronic  diatancea  in 
the  Integrand. 


3. 


II,  ?OBMUUTION  OF  THX  PROBLEM 

Without  using  carsfully  deslgnsd  foms  of  the  correlated 
wave  functions  we  would  have  to  evaluate  Integrals  of  a  conald- 
arabla  coaplaxlty.  But  If  we  e.g.  use  trial  wave  functions  of  the 
fonn 

where  the  analytical  fora  of  ^  and  ^  are  restricted  to  certain 
laportant  als^tla  cases.  It  Is  possible  to  aake  progress  in  the 
calculations  of  the  occurring  Integrals,  ^  aay  be  an  anttsynae- 
trlaad  product  of  suitably  chosen  spin  orbitals  or  a  superposition 
of  such  products,  aiid  ^  a  nodsless  function,  syanetrlcal  with 

,  r  u  ,  ,  a, 

respect  to  perautatlons  of  the  electronic  coordinates. 


9}  L.  Redel  and  P.0,  Ljh^ln,  Phys,  Rev.  114  ,  752  (  1952} 


Is  here  a  coablned  space-spin  variable,  such  that 

Is  the  position  vector  of  the  electron  \  with  respect  to  the  nucleus 

and  a  spin  variable.  Slightly  aore  general  trial  functions  of  the 

saas  type  have  also  been  suggested  with  ^  not  necessarily  a  function 

31 

of  all  the  Interelectronlc  distances. 

where  and,^  Is  the  antlsyaaetriser  in  M-partlcle  space. 

The  correlation  factor  ^  has  often  been  chosen  as  a  polynoalal  of 
Y7  and with  suitable  adjustable  paraaaters,  but  also  other 
slopla  analytical  foras  have  been  used.^^’^^^ 


4. 


10)  P.  Ualah  and  S.  Borotfitz,  Phys.  Rev.  119  ,  1274  (I960} 


If  we  keep  In  mind  all  these  end  other  cooBion  restrictions 
on  the  correlated  wave  functions  used  so  far  in  atomic  calculation 
it  may  be  general  enough  to  start  out  by  treating  an  integral  of  the 
form 

This  may  factorize  into  integrals  of  lower  dimensions  where  the 
Integration  variables  are  coupled  through  the  Interelectronic  dis¬ 
tances.  The  simplest  non-trivial  form  of  such  a  factor  would  be 


where  we  have  used  the  orbital  form  explicitly. 

Here  *  spherical  harmonic  as  defined  e.g.  in  Bdsunds'  book. 


11)  A.R.  Bdsunds,  "Angular  Homentum  in  Quantum  Mechanics",  (Princeton 
University  Press,  1957) 


If  we  now  expand  a  product  of  spherical  harmonics  according 


to  the  formula 


(5) 


{-Ip  p 

wh«r«  vn. «  )lt, -  Vtt^  ,  we  can  write  i\  aa  a  aua  of  Intagrala 

-  %  { 

•Ut)  J^,t';'{vi>;)f.(>;) 


(6) 


For  inatance  K(aoo,  410,  cio|<ioo,e.i0^^lo)  niii 

T«*0  ■^•OP  -^PPP  ^OOP 

•0» ,  1.  .  „  1  .  «n<l  1  .  . 

^  II  ^  f  OXO>  ♦OX 

Hare  Ilf;)  ^l'9)s 

and  ^  ^  TTe  tn-^-vn^  . 

Tbe  Higner  3-J  ayad>ol8  are  defined  in  reference  11.  Further  we  can 
wrltal2) 


12)  F.O.  Lowdln,  "A  Theoratlcal  Inveatlgation  Into  Soae  Fropartlaa 
of  Ionic  Cryatala**,  Theala,  Dppaala.  1948  (Alwiviat  and 
Hikaall) . 

-f*  Coaipara  fonaila  (4.6.5)  In  raferance  11,  and  apeclflc 


foraulaa  In  reference  15. 


6. 


IV 


^A*f  ^:*^;  -  Zl^K  C04  l)^j  j  1*  th«  Lvgvndr* 


‘1  ^  ■  I  ~ 

polynoalal  of  tfogroe  tV  and 


Bquatlon  (7)  Cogathar  with  tho  addition  thooroa, 

(9)  (coifif  )  s  Alt  ( Z  »1.+  y)'iL  yJiitvy.jpi) 

«  '  VH—IV  '  ' 


than  givo 


(10) 

whom 

(lOa) 


and 


4iu.'v  __ 

’  -W,i,nv  ^ 

iMu^n  0} 

^  90  |4 


7. 


Vfc‘ 


(10b) 


n  K  n 

j.--li  r.-l  -t.  -itw 

ft  zrt.  ^ 

§  A 

*  V'  # 

By  chmga  of  intogratlon  varlabl*  ,  where 


,  we  get  froai  equation  (8) 


(11) 


*  !ff«(j|  *■ 

Depending  on  the  analytical  font  of  the  functlona  ( 1^)  i 

we  can  now  proceed  In  different  ways.  We  can  e.g.  calculate  the 
^l^i)  uelng  the  ’'Q-polynoaiala"  of  Lowdln.^^^Thla 


13)  P.O.  Lwdln,  Adv.  In  Phya.  S  ,1  (1956);  ’'Quantiaa  Theory  of 
Cohesive  Properties  of  Solids."  p.98 


aisans  that  we  can  use  the  well-known  explicit  expression  for  the 
Legendre  pelynoodals  and  write 


(12) 


— .  .  K*K 

“tiji't.’llji) *  (2^c* 

Here  Q>  (.^.1^)  is  a  hoaogeneous  polynoid.al  In  k'^  and  ^ 


8. 


of  dograa  .  Sacurslon  fonmilaa  for  tha  are  given  in 

rafaranea  13.  Tha  axpraaaions  wa  gat  for  in  this 

way  my  ba  tadloua  and  lengthy  to  work  with  but  are  in  principle 
not  co^licatad. 

TO  ba  apaelfic  wa  have  chosen  to  treat,  in  the  following, the 
casa  iriian  ^  la  agtprasslons  (1)  and  (2)  is  a  polynomial  in  K  , 
and  .  Without  loas  of  generality  wa  can  than  choose  the  functions 


in  tha  Integrals  (3)  and  (4)  to  simply  be  the  powers  of  the 
intaralaetronle  distances.  From  equations  (11)  and  (10a)  we  get 

t  '  Xvl  t  5 H-+I 


(  *  _  V>«-1 

It  is  now  in  general  necessary 


where  X«  •  * 

to  distinguiah  between  tha  regions  and  This  splits  U^^^nto 

six  integrals  with  different  order  relations  between  the  integration 
variables,  i.a. 


f  \  " 

•  a  1  ' 

**  7  '* 

••  «  •* 


9. 


M  M 


j.  J.l'"'].  i  j  j '  ^ 


H«r« 


(15.) 

(15.)  t:41«>i)=-vac 

i-o  a 

«h«r«  [  ^] 

(15c)  ®  0  for  -t  i  m  +  I 


<♦  I 

*  Integer  part  of  ,  and 


(15.) 

Equation  (15c)  showa  that  if  one  of  X  ,  p,  and  Is  an  even 
nuaibar,  one  of  the  three  infinite  sianatlons  in  equation  (10)  is 
reduced  to  a  sun  over  a  finite  number  of  terms. 

In  order  to  proceed  further  with  calculation  of 
we  have  to  be  more  specific  about  the  form  of  the  radial  part  of 
the  orbitals.  We  are  in  the  following  going  to  limit  the  discussion 
to  orbitels  of  the  form 
(16)  twir).  ejcr^>(-av  ■Jr'' 

where  is  an  integer.  This  allows  us  to  include  in  our  treat¬ 

ment  radial  parts  expressed  as  hydrogen- like  functions,  Laguerre 
functions^^and  linear  combinations  of  Slater  or  Dirlcblet  functions! 


14)  0.0.  Hall,  Rap.  on  Progr.  in  Phys.  Vol  XXII  (1959)  p.ll 


10. 


In  then*  €••••  th«  from  aquation  (14)  that  we  want 

to  calculate  would  be  auparpoaltlona  of  the  ’**  8*^ 

from  ualng  tha  apactflc  form  In  aquation  (16)  aa  radial  parts. 

Wa  now  turn  to  the  problem  of  calculating  in  coefficients 
defined  in  equation  (10b) .  Here  we  have  to  in¬ 
tegrate  over  products  of  three  spherical  harmonics.  In  order  to 
get  a  non-vanishing  contribution  from  an  integral 

it  is  well  known  that  the  following  relations  have  to  be  fulfilled 

■«“* 

These  restrictions  give  among  other  things  that 


(17) 


‘J  Y  u)y(2)  y(2) 


and  that  (,-^4  j  Vm.  |  ij-r )  *  0 

or  -tw+v dZ*!.  ^  ..  .  )  Schweinler^^^has  given 

useful  explicit  formulas  for  expanding  products  of  spherical  harmonics 
of  the  sasw  arguawnt  in  spherical  harmonics.  Formulas  of  that  typo  are 

feaaable  to  use  in  calculation  of  the  integrals  in  (17). 


15)  H.C.  Schweindler,  Quarterly  Progress  Report,  M.I.T. ,  Solid 
State  and  Molecular  Theory  Group,  January  15  (1954)  p.51 
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III.  FORMULAS  F(»  CALCULATION 


In  the  nuaerical  evaluation  we  have  extensively  used  a 
set  of  three  auxiliary  functions^^defined  as 

to 

(•)  e.* ’■  cU- 

••  *• 

(18)  (b)  j 

pO  00 

(c)  A)* j  I  i  ^  ^ 


We  will  first  discuss  the  “spherically  synnetric  case"  ^ „ 

This  we  hope  will  show  the  owin  idea  of  the  method  and  the  general¬ 
ization  to  other  cases  is  straightforward.  in  equation 
(10b)  will  then  reduce  to 


(19)  JB^^^Uoloojoo)  «  ')• 

This  together  with  the  restriction  on  the  form  of  orbitals  used,  as 
expressed  in  equation  (16)  gives 

to 

I  .  -ax  ^ 

where  defined  as  in  equation  (13)  with  x' 

*”**  ^5(a)=C.^  i*  .  It  should  be  noted  that  we 

have  not  normalized  the  orbitals. 

The  integrals  I  K  1  ^  *51  ^  1  ^  p  )  can  now  be  expressed 

in  the  auxiliary  functions  A  ,  N/  and  \j(/  .  We  have  e.g. 


12. 


^■‘Ixoo  {>+ii 

*H  (xJ,)  +V(^»c|<|U*2i^^*>*2.2v)J 

^**  r^i]  [Kii] 

I(*/»)f|^4*|>(.o).i(J..i)((..»)>l£l  C  (aV.‘,)(ll'f) 

i^A  |iA  * 

(21b)  ♦*+X  ‘►^•■14  -Z^  \+  W(*(|^Y‘*t«2i^2|j  ii) 

•od  M, 

V'  > 2-4) 

+  y/(^][Kl^»X<2-lk^  b+1,  ^*4)  ••■  '<V 

Both  (21«)  and  (21b)  rapraaant  may  Integrals  of  interest  in 
atosdc  theory.  In  particular  (21a)  for  X"  land  A(||U^splaced  by  I 
can  be  used  for  calculating  the  ordinary  two-electron  repulsion  Inte¬ 
grals  for  spherically  sysnetrlc  functions.  The  Integral  In  aquation 
(22)  Is  tha  slaiplast  of  the  1  I  X|4  v)  U  that  requires  an 

Infinite  sum  of  .  in  connection  with  (21a)  and  (21b)  we  should 


■Lao  note  the  trivial  but  uaeful  relationa  like 


For  certain  restrictlona  on  the  fom  of  the  correlated  wave 

function  uaed|Hhe  aost  coapllcated  Integrals  bound  to  occur  would  be 

_  4mv 

possible  to  express  as  linear  conbinatlons  of  J 


16)  1.  Ssasx,  Quarterly  Progress  Report  M.X.T.  Solid  State  and 
Molecular  Theory  Group,  July  15,  1962 


To  show  how  such  an  Integral  can  be  calculated  we  are  In 
section  V  and  In  the  Appendix  discussing  In  particular 


*9 

* y  f ‘t .  1 .  )  ■►w  y  a  ,p )  ♦  w  ( y  I  4 


♦ 


U5) 


t  Wk)4V(|kKy\^4*Jk,^0>-2Vo) 

♦  W  (y^  ic,  ,  j,-2  -i. » ,^-2Vc) 

I  ,s,*2-2to) 

v*2*2to,^*5,p*2-2W)tW(ejr^Jj»4i»»2»c,  s*5,  ^-2tej 

y  ^  *i'*'  I  * .  p*i'"2-k)*W(ye^|  i*2*2k, ,  ^5,  k)\ 


14. 


“  (24cO)  (Zk-l) '  v2ic)-vW(yiOtyl<^«Z*2<c,f*5,  vik) 


+  W(y^fli\*4M42k,<^»l,^*2-2k)*W(*(y^lp4i<42V<,s«\,<y*2  >ik) 

*  W{^^e(  j<^*i»2lc.  ■f>«s,p*2-2lcj*vi/(^C(^is*i»-*2»Cj  p»  l,^*2-2  tc)^ 

*W  <^*2‘2k,  &»  I,  f  k-2»c)*W(jf*|ilw-2lc.,p*l,^»2-2k,)j^ 


All  these  formulas  do  not  tell  very  much  about  the  method 
unless  we  also  describe  how  to  calculate  the  A4>  V'^  ,  and  W/i 
and  Investigate  the  convergence  of  the  Infinite  sums.  This  will  be 
done  in  sections  IV  and  V. 

If  at  least  ot^o{X\  ,  \  and  ^  in  equation  (10b)  is  different 
from  zero  we  will  get  more  complex  situations.  For  example  , 

Ul  giveB^.^\»Ojlol(0H.l  and3^^(^«0''0|'0)=-t< 

as  the  only  nonvanishing  coefficients.  It  is,  Wowever,  also  for  a  gen- 


—  4».  w 

eral  ,  a  straight  forward  procedure  to  express  it  in  terms  of 

A'a  ,  Vv^  And  W'/),  by  going  back  to  formula  (14).  This  means  that  al¬ 
though  we  in  sections  IV  and  V  are  practically  treating  the  spherically 

IOtf  O 

>  i»B*=hod8  and  results  from  there  are  with  very 
few  exceptions  directly  applicable  to  the  general  case 


15. 


IV.  COKPtmVTlONAL  SCHEME 


In  calculating  a  particular 
have  to  dlstlngulah  between  two  casea,  -tt^O  andii^O. 

1.  When 'ttiO  the  calculation  la  aet  up  In  three  atepa. 

(1)  First  we  fora  a  table  of  all  for 

j t  using  the  relation 


we 


(24) 


where  i  I  . 

(11)  Frou  this  table  we  can  then  construct  a  table  of 

using  the  foroula 

(25)  Kl^,-i!c.^«ys-ru)  { AC«y^+j[ 

which  la  valid  for  4^0  and 


(Ill)  The  flnel  step  la  now  to  use  the  VI4  to  calculate 
the  W(  I  Interest.  This  can  be  done  with  the 


recurrence  relation 


(26)  + 

for  froa  0  4o  4w  .  (26)  Is  not  valid  unleas^^O 
+  i  - 1  . 


16. 


1^99 

0,0  that  ara  expreated  as  infinite  sums  like 
those  in  equations  (22)  and  (23)  we  need  in  addition  to  calculate 
^  ^ -ti}  for -U  40.  This  can  not  conveniently  be  done  with 
(24) -(26),  partly  because  it  is  a  tedious  procedure  to  find  all  the 
“starting  valuas“  I  )  and  -l}  that  are 

necessary,  and  partly  due  to  the  fact  that  these  recurrence  relatione 
will  in  this  ease  not  be  numerically  stable,  i.e.  loss  of  significant 
figures  is  bound  to  occur.  In  this  case  (4t4  6)  we  use  the  formula 


(27) 


idiere  for 


M 


.H-l..  •  •  I 


and  for 


and  nJ  are  in  principle  Infinite,  but  in  ectual  calculations 
they  are  of  course  assigned  finite  values.  In  order  to  avoid  very 
large  and  very  amsll  quantities  these  forsailas  have  been  combined  to 
the  following  coeputational  schesm: 


(30) 

where 

and 


17. 


^  .„d  5,  .  1 

(^*1 )  ('l-^^^H  +  V)  ' 


U*a) 


«d  X.. 


.,lwlth  M-a 


For  calculations  of  ^  *  computational 

schaaM  may  be  darlvad  from  (28) 


(31)  Z '^y 

^here  ^..  with  'VT; 

*  (p<t^y)  (y*‘)  +  m+v>)  '  • 


,(K*A+y  -t+i+A+H 


18. 


V.  CONVBRGBMCB  PRORRTXE8 

Hi  Th«  achana  for  calculation,  of  W  with  'tt<0  , 
defined  in  (30)  Includea  three  infinite  auBM  and  the  rate  of  conver¬ 
gence  for  each  of  th«B  la  eaaantial. 

In  actual  calculatlona  t(,  and  ^  will  be  of  roughly  the 
same  order  of  augnitude  and  we  will  lladt  our  diacuaaiona  of  the 
convergence  propertlea  to  thia  caae. 

The  factora  containing  |it  or  ^  in  the  recuralon  foraulaa  for 
<|^and  ,  converge  to  unity  and  the  latter  nay  be  approximated 
aa  followa 


(32) 


i.* 


Obvioualy  the  convergence  for  )/^and  ^ia  faater  than  for 
^  in  (30)  la  the  critical  part. 

Actually  ^la  zero  for  H  infinite,  and  (30)  may  be  replaced 


by 

(33) 


M 


The  nuad>er  of  multipllcationa  (or  dlviaiona)  needed  to  calcu¬ 
late  W  by  (30)  and  by  (33)  ia  reapectively 


and  •  tM  , 


(34) 


lOH+ihi 


19. 


R«elMr  th«n  giving  «n  extnnslva  thcorncical  traatawnt  6t  the 
convergence  properties  of  W  ,  we  have  chosen  to  exeag>ll£y  its 
behavior  by  calculating  W  |  b  0  -  I  )  for  some  typical  values 
for  .  As  seen  from  formulas  (27)  and  (31)  the  choice  of 

^  ~  )  represents  some  of  the  worst  convergence  prop* 

erties  possible.  In  table  I  wa  have  tabulated  the  nunber  of  signi¬ 
ficant  decimal  figures  dt  of  W  .  calculated  by  (30),  as  a  function 
of  H  and  hi  .  Table  II  shows  the  number  of  significant  figures  d* 
as  a  function  of  M  ,  when  W  is  calculated  by  (33) . 

From  (34)  and  Tables  1  end  II  algorithms  may  be  derived  for 
the  choice  of  M  and  nI  .  In  principle  one  could  even  give  M  and  hJ 
as  explicit  functions  of  the  arguments  desired 

number  of  significant  figures  d. 


[  I  I  >1^  As  an  example  of  the  behavior  ofW  with 
negative  ,  we  have  chosen  to  treat  X,,,  for  V* 

The  explicit  fonouAa  for  calculation  of  the  Integral  is  given  in 
(23),  but  to  simplify  the  discussion  ws  rewrite  this  eq^uatlon  as 

05)  l-l  )  'Xt 

where  K  la  In  principle  infinite,  end  the  definitions  of  J)^andl)^ 
are  obvious. 

For  each  of  the  fourteen ^^'tappearing  in  3^^  ,  we  have 


the  relation 


Jt4v»  *.  +  7 


20. 


Whcn*^  Is  Isrgs  enough  all  the  ^*4  will  have  -^<0  end 
W,  in  (30)  %rlll  be  approximately  constent.  Further 
la  obviously  independent  of  -k  ,  so  even  under  the  summation  over 
the  variation  in will  be  represented  by  the  factor  (-^4)) 
and  this  variation  is  very  small  for  large  .  In  fact  this  is 
the  main  variation  for  each  of  the  .  The  approxlmste  constancy 

of  the  W'A  ,  shows  that  the'^^^  behave  like  'X. 
a  good  estimation  of  the  convergence  properties  of 


and  this  gives 

•• 

-<(r| 


21. 


Table  I 

The  number  of  eignlf leant  tiguree, 
calculated  by  (30). 
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We  have  chosen  to  carry  16  decimal  figures  (Fortran  double  precision) 
In  our  calculations,  because  in  combining  the  X‘A  to  matrix  elements, 
necessary  in  the  variational  calculation,  cancellation  of  significant 
figures  will  occur  In  some  cases  (e.g.  when  the  radial  parts  of  the 
orbitals  are  Lagaarre  functions) . 


22 


T«bla  II 

. .  “  I 

Th«  nuaber  of  significant  figures  In  )  0,0^- l} 

calculstad  by  (33). 


200 


14 


14 


23 


mEMDIX 


At  tn  MUaplt  w*  will  dltcutt  tht  Inttgrtlt  produced  by  a 
correlated  wave  function  for  a  three  electron  atonic  problem.  An 
adequate  trial  wave  function  for  a  doublet  etate  can  be  written 


(Al) 

where 

*ni 


After  Integration  over  the  spin  varlablet  the  average  value 
of  the  Haniltonlan  bdconet 


where  it  the  permutation  operation  that  interchanget  the  electron 
coordinate!  2  and  3.  For  an  S-atate  we  nay  chooae  the  apace  functiona 


(A2) 

and 


•§^13)  »  >s 


The  only  part  of  tha  Haad.ltonian,  that  Involvea  calculation 
of  an  integral  with  \  ^  and  V  all  0 


ia 


The  evaluation  of 


24. 


^ +2cX(2,A.^2,A»^yi4|6  0Zl  l-l  0 

+c*i(2.a.,2<l,2ilOOlU-10)'lUt4,^<t+4,Z«i|llo|-IDO) 

vhare 

«  1  (2ft., z«., 24|o02.jo00)  +  iclUft., 2.^1, z4loozlioo) 

+  cl(2ft,Z(;,,^4lo02|200)-l(24,AU,<l+4loillooo) 
-2clC2A.,a.*4,<i^tloii|  loo)-  c*lC2^<^a.+i,ft4|DI||ilc/) 


will  serve  to  deiDonatrate  the  convergence  of  the  Infinite  series  In 
(23)  for  celculetlon  of  ^  I  ^  I  *  ^  3  elso  the  lapor- 

tance  of  this  term  In  the  average  value.  For  the  Ll-atom  <t*2.tC  , 
and  C'O.'hb  would  give  a  fair  approximation^^  to  the  ground 
state  wave  function. 

In  Table  III  we  list'b,  end  the  first  fifteen  In  the 
Infinite  sum,  that  appears  in  the  definition  (35)  and  (23)  for 

l(2a  tc*  ^loU  ht'*' I  J  •  observed  convergence  for  the 

'  '  *'*  -J, 

series  Is  actually  better  than  for  ^  ^  and  It  Is  Justified  to 

■k«i 


conclude  that  the  reminder  in  general  is 

*4 


E  A  ^  ^  A 


25. 


ifhleh  la  our  mum^Io  aaant  that  «•  bava  aisht  algnlf leant  flguraa  In 

X. 

Iha^A  appaarlng  In  aaeh^^^  vara  aalculatad  by  (30)  irlth 
M*l4a4o,  uhleh  aeeordlng  to  Tabla  X  glvaa  alavan  or  twalva  algnl- 
floant  flfuraa*  aai  aa  no  appraelabla  eaneallatlons  occur  In  tha 
foraatlon  of  aaehl^thla  doaa  not  Influaaca  tha  accuracy  of 

Tabla  If*  llata  tha  IntagralaX  and  «a  Ilka  to  point  out 
tha  tha  aboua  aantlonad  Intagral  l)l/'l)tnay  da- 

flaltaly  not  bo  naglaetad  In  tha  calchlatlon  of  tha  avaraga  valua. 
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Table  HI 


for  k  s  0,  1,  Z,  ...  12,  as  defined  by  (23)  and  (35)  for  the 
integral  1(5.72,  4.26,  4.26  J  0  1  1  |  1  1  <l) 


Dq  »  0.358  027  6958  •  10"® 

Dj  o  0.004  913  4950  •  lO"® 

D2  *  0.000  046  7536  •  10"® 

D3  a  0.000  002  9199  *  10*® 

=  0.000  000  3819  *  10"* 
Dg  *  0.000  000  0759  •  10*® 

=  0.000  000  0199  •  10*® 
=  0,000  000  0063  •  10*® 
Dg  =  0.000  000  0023  •  10*® 

=  0.000  000  0010  •  10‘® 
Djq  =  0.000  000  0004  •  10’® 

Djj  =  0.000  000  0002  •  10*® 

Dj2  =  0*000  000  0001  •  10*® 


and  consequently 

1(5.72,  4.26,  4.26  1  0  I  1  |  I  1  - 1 )  =  0.  362  991  35  •  lO'®  . 


Errata 


Formulas  (23).and  (35)  for  calculation  of  Ip9*I  1  ^  ~1) 

are  not  correct.  The  term 


0, 

should  be  deleted  and  the  summation  ^  replaced  by  iZl  •  The 
correct  formulas  were  used  for  the  cafculattoni  reported  l^^Tables  111 
aud  IV. 


Table  IV 


Integrals  I(a,  p,r|p»q,  s|Xt|i>  v)  appearing  in  formula  (A3). 
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A  =  1.930  928  02  •  10'® 


0. 587  090  36 
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